Superfiuid-insulator transition of disordered bosons in one-dimension. 



0^ 
O 
O 
(N 

(N 

C/3 



Ehud Altman^, Yariv Kafri^, Anatoli Polkovnikov''. Gil Refael"* 
^Department of Condensed Matter Physics, The Weizmann Institute of Science Rehovot, 76100, Israel 
"^Department of Physics, Technion, Haifa 32000, Israel 
^Department of Physics, Boston University, Boston, MA 02215 
'^Dept. of Physics, California Institute of Technology, MC 114-36, Pasadena, CA 91125 

We study the superfluid-insulator transition in a one dimensional system of interacting bosons, 
modeled as a disordered Josephson array, using a strong randomness real space renormalization 
group technique. Unlike perturbative methods, this approach does not suffer from run-away flows 
and allows us to study the complete phase diagram. We show that the superfluid insulator transition 
is always Kosterlitz- Thouless like in the way that length and time scales diverge at the critical point. 
Interestingly however, we find that the transition at strong disorder occurs at a non universal value 
of the Luttinger parameter, which depends on the disorder strength. This result places the transition 
in a universality class different from the weak disorder transition first analyzed by Giamarchi and 
Schulz [Europhys. Lett. 3, 1287 (1987)]. While the details of the disorder potential are unimportant 
at the critical point, the type of disorder does infiuence the properties of the insulating phases. We 
find three classes of insulators which arise for different classes of disorder potential. For disorder only 
in the charging energies and Josephson coupling constants, at integer filling we find an incompressible 
but gapless Mott glass phase. If both integer and half integer filling factors are allowed then the 
corresponding phase is a random singlet insulator, which has a divergent compressibility. Finally in 
a generic disorder potential the insulator is a Bose glass with a finite compressibility. 
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I. INTRODUCTION 

Superfluid-insulator transitions occur in a variety of 
experimental systems, ranging from low-temperature He- 
lium through Josephson arrays to ultra-cold atomic sys- 
tems. The simplest paradigm of such a transition is 
the rather well understood Mott transition of interact- 
ing bosons on a perfect lattice commensurate with the 
boson density. [1, 2]. The theoretical picture is far less 
clear in disordered systems, which occur in a wide variety 
of experiments: Helium in Vycor, superconductor-metal 
and superconductor-insulator transitions in nanowires 
and thin films [3-15]. Recently, disordered systems were 
also realized using ultracold atoms [16-21]. Furthermore, 
this topic was brought back into the limelight with recent 
experiments in solid Helium-4, which show the appear- 
ance of a superfluid fraction [22] . One suggested expla- 
nation for this phenomenon is Helium turning superfluid 
in structural defects of the surrounding solid [23] . 

Of particular interest is the superfluid-insulator transi- 
tion in disordered one-dimensional systems. Even with- 
out the disorder the superfluid phase in one dimension 
is more subtle than in high dimensions. In particular it 
does not exhibit true long range order. Nevertheless the 
uniform superfluid admits a simple description in terms 
of a universal harmonic theory, or Luttinger liquid. In 
the opposite limit of a disordered potential but no inter- 
actions, particles are always localized. One might naively 
guess that there is no superfluid phase in the presence of 
disorder since interaction alone or disorder alone both 
have a localizing effect on the bosons. This however does 
not seem to be the case. 

The simplest way to see this is to introduce disorder 
as a perturbation to the interacting superfluid within the 



Luttinger liquid description. This was done by Giamarchi 
and Schulz in Refs. [24, 25]. The main result of this 
approach is to describe a phase transition between an 
essentially uniform superfluid, in which the disorder is 
irrelevant, into a localized phase. The natural tuning pa- 
rameter of the transition is the interaction constant and 
it occurs at a universal value of the Luttinger parameter, 
independent of the strength of the disorder. 

The above approach suffers from two main limitations. 
First, because it is perturbative in the disorder strength 
localization is signaled by a runaway RG flow. There- 
fore the approach does not allow for a detailed theory 
of the insulating phase. Second, the natural regime for 
the phase transition in this analysis is that of strong in- 
teractions and a nearly uniform superfluid, which is not 
always the case in systems of interest. For example atom- 
chip traps, in which ultracold atoms seem to undergo a 
localization transition [17], are in precisely the opposite 
regime. The bosons are weakly interacting, while the po- 
tential they feel is highly disordered [26]. It is not clear 
whether the analysis of Giamarchi and Schulz provides a 
valid description of the transition in such a system. 

Different approaches have been used to specifically de- 
scribe the insulating phases of bosons and suggested sev- 
eral possibilities depending on the nature of the system. 
In the most generic disordered potential, Ref. [27] argued 
for the formation of a Bose glass phase characterized by a 
finite compressibility and diverging local superfluid sus- 
ceptibility. In the presence of a commensurate lattice 
Refs. [28, 29] predicted the existence of a Mott-glass 
phase, an incompressible yet gapless insulator. 

In recent work we introduced a unified approach to 
treat both the phase transition at strong disorder as well 
as the properties of the insulating phases [30, 31] . For this 
purpose we employed a real space renormalization group 
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(RSRG) technique [32-35]. We found that the super- 
fluid insulator transition at strong disorder is insensitive 
to the type of disorder introduced into the system. It 
is always Kosterhtz-Thouless like in the following sense: 
characteristic time scales and length scales both diverge 
at the transition as exp(l / — ac), where a is the tun- 
ing parameter. The nature of the disordered superfluid 
phase is also universal. It is described by an effective 
harmonic chain with random Josephson couplings drawn 
from universal distributions generated as fixed points of 
the RSRG flow. These distributions were recently used to 
compute the localization behavior of density waves [36]. 

The symmetry properties of the disorder, while not im- 
portant in the superfluid phase or the transition, are cru- 
cial for determining the nature of the insulating phases. 
Using the RSRG approach we confirmed the formation 
of a Bose glass phase for generic disorder and a Mott 
glass for a commensurate lattice with off-diagonal dis- 
order. The latter phase was also seen in recent numer- 
ical simulations [37, 38]. In addition we found a novel 
glassy phase, which we termed a random-singlet glass, 
in a system with particle hole symmetry. This phase 
is characterized by a divergence of both compressibility 
and superfluid susceptibility. Nevertheless it is still in- 
sulating, with conductance dropping as exp(— \/L) with 
length. This phase is analogous to the random-singlet 
phase found in the spin-^ X-Y chain [33]. 

The purpose of the present paper is twofold. First, we 
provide the detailed analysis of chains with generic di- 
agonal disorder, leading to the results of Ref. [31]. Sec- 
ond we extend the analysis and compute the value of 
the Luttinger parameter at the phase transition within 
the RSRG method. We find that, at strong disorder, the 
transition occurs at a non universal value of the Luttinger 
parameter that depends of the strength of disorder. This 
is contrary to the perturbative analysis of Refs. [24, 25]. 

The structure of the paper is as follows. In section II 
we define the model we study and discuss its relevance to 
actual physical systems. We give a detailed derivation of 
the RSRG flow equations for the special case of particle- 
hole symmetric disorder in section III and for generic 
disorder in section IV. We give a detailed account of the 
numerical as well as the approximate analytical solutions 
of the flow equations. Then in section V we solve for 
the value of the Luttinger parameter at the transition. 
Finally we conclude with a summary of the results and a 
discussion of their possible experimental implications. 

II. MODEL 

Our starting point for the theoretical analysis is the 
quantum rotor Hamiltonian 

n^^Vj {fij - Hjf - Jj cos {ipj+i - ipj) . (1) 
i i 

This model describes an effective Josephson junction ar- 
ray with random Josephson coupling Ji and charging en- 



ergies Ui. In addition there is a random offset charge 
to each grain, which is tantamount to a random gate volt- 
age. Although the model can be visualized as a Joseph- 
son junction array, it actually provides an effective de- 
scription valid for a wide variety of systems that undergo 
a superfluid to insulator transition. In bosonic systems in 
particular, such transitions arc usually driven by quan- 
tum phase fluctuations. The hamiltonian (1) should then 
be thought of as a low energy effective theory one obtains 
after integrating out the gapped amplitude fluctuations. 
The remaining degrees of freedom relevant to the transi- 
tion are the quantum rotors. 

One concrete example of how such a model is natu- 
rally generated at low energies is provided by a system 
of ultracold atoms in an atom-chip trap. In this system 
the disordered potential is induced by corrugation in the 
wire that generates the trapping magnetic field[26]. With 
increasing corrugation, the atoms concentrate in small 
puddles at minima of the potential. Neighboring puddles 
are connected with each other by a random Josephson 
coupling which depends on the potential barrier between 
them. The result is exactly the random Josephson array 
defined in Eq. (1). 

Another possible physical realization of the model 
(1) is a disordered superconducting nano-wire. Here 
the issue is more subtle because there may be gapless 
Fermionic degrees of freedom that generate dissipation. 
Indeed Refs. [39, 40], apphed the RSRG to such wires 
starting from a Hertz-Millis[41, 42] dissipative action, 
with a dissipation term |a;|'0*^/;. An alternative approach 
is to use phase only models, which describe resistively 
shunted Josephson junction arrays[43]. This naturally 
leads to a dissipation term of the form q'^\u!\(j)* (j), which 
docs not affect global superconductor insulator transi- 
tions. Such models combined with strong disorder may 
also be described by the present analysis in parts of their 
phase diagram. 



III. PARTICLE-HOLE SYMMETRIC 
CHEMICAL POTENTIAL DISORDER 

Of all the random coupling constants in the model (1), 
the random offset charge (or local chemical potential) 
seems to be the hardest to incorporate in an RG treat- 
ment. Since the offsets simply add up to give the offset of 
a block of sites, it seems clear that this disorder will just 
grow as the square root of the scale of the real space RG 
making it hard to track. However this difficulty turns 
out to be largely superficial and can be easily overcome 
in the analysis. To make the discussion more transparent 
we start from the case where only integer and half-integer 
offset charges, rii, are allowed on each site. This con- 
dition maintains particle-hole symmetry, and therefore 
still does not correspond to the generic case. Neverthe- 
less, this restriction allows a relatively simple RG analysis 
which affords important analytic and numerical insights 
into the possible phases and the phase transitions. In the 
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next section we generalize our treatment to the case of 
generic disorder. 

We note that despite the restrictive condition, allow- 
ing only rij ~ 0, 1/2, this type of disorder may actually 
be a reasonable approximation for chains of supercon- 
ducting grains with pairing gap much larger than the 
charging energy. Under these conditions we can assume 
that the electrons on the grains are always paired and 
we can take e* = 2e as the unit of the bosonic charges. 
On the other hand the positive background charge, is a 
random number that could be even or odd in units of e 
and consequently either integer or half integer in units 
of the boson charge e*. Allowing for charged impurities 
on the substrate or unscreened coulomb interactions be- 
tween different grains, would of course lead violation of 
the restrictions on the off set charges. 



A. Particle-hole symmetric quantum rotor model 

The essence of the renormalization group transforma- 
tions cither in real or momentum space is the gradual 
coarse graining of the system. In this section we extend 
the decimation scheme of Ref. [30] to the Hamiltonian in 
Eq. (1) for the case that rij can take the values of and 
1/2 randomly. The last condition ensures the particle- 
hole symmetry in the problem: the Hamiltonian does not 
change under the transformation n ^ 1 — n. These two 
values of Uj represent the two possible extremes which 
drive the physics of the Bosc-Glass [27]. Sites with rij = 
have a well defined Coulomb blockade with charging en- 
ergy Uj. Sites with rij = 1/2, on the other hand, have no 
Coulomb blockade. With no further interactions, these 
sites yield both infinite compressibility and infinite super- 
fluid susceptibility due to the number fluctuations costing 
no energy. The Hamiltonian (1) is characterized by the 
distribution of hoppings Jj and charging energies Uj , and 
of the proportion of sites with n = 1/2. We will refer to 
the latter sites as half-integer sites or 'half-sites'. 

B. Extended real-space renormalization group 

Let us now construct the extended decimation scheme 
for the model (1). Following Refs. [30, 32-34] we con- 
struct an RG scheme that eliminates iteratively large en- 
ergy scales from the Hamiltonian. Two sites connected 
by the strongest bond will be converted to a phase- 
coherent cluster. Similarly, in sites with strong charging 
energy U we eliminate all the excited states. However, 
the result of this elimination will be different for integer 
and half-integer sites. Let us now discuss these steps in 
detail. 

We denote the largest energy scale in the Hamiltonian 
(1) fi = max{ J,;, J7i}. In each step in the RG we elim- 
inate the strongest coupling from the Hamiltonian, and 
hence successively reduce fl. If the strongest coupling 
is the charging energy of site i, Ui, we eliminate all the 



excited states of this site. For integer sites with n.; = 0, 
we minimize the charging energy by setting rii — and 
include the coupling of this site to the rest of the chain 
perturbatively. As in Ref. [30], the second order pertur- 
bation theory leads to a new coupling between the new 
nearest neighbors i — 1 and i + 1: 

Ji-l,i+l = Ji-lJi/^- (2) 

On the other hand, if rii = 1/2 we reduce site's i Hilbert 
space to the states = 0, 1. The hoppings connecting 
site i to its neighbors are still active, and to the first 
approximation are not affected by the elimination of the 
high energy states. The decimation step for rii = 1/2 
produces a new kind of site, a doublet site, only capable 
of having = or 1. Let us denote the fraction of 
doublet sites as s, the fraction of integer sites as q, and 
the fraction of half-sites as p. Note that these fractions 
add up to unity p -\- q + s = \. 

When the strongest coupling in the chain is the bond 
Ji, unless both sites i and are already-decimated dou- 
blet sites, a phase-coherent cluster forms. Since charging 
energy is the inverse of capacitance, the effective Ui^ i+i 
of the new cluster will be: 

1 _ 1 1 

c^^ caTT ^ ' 

For a doublet site, Ui is set to c». It is easy to see that 
the filling factor n is an additive quantity: 

n.i^i+i = {ui +rT,:+i) mod 1. (4) 

Therefore two half-sites or two integer-sites form an 
integer-cluster. An integer site and a half-site form a 
half-cluster. Similarly, a doublet site and an integer site 
form a half-cluster, and a doublet and half-site form an 
integer cluster. 

It is important to note here that the above decima- 
tion step does not assume long range order; it states that 
phase fiuctuations within the newly-formed cluster are 
harmonic, and therefore the cluster can not be broken 
due to phase-slips. These harmonic fiuctuations are cru- 
cial for the understanding of the properties of the super- 
fluid phase, as explained in Sec. V. Nevertheless these 
phase fluctuations can be neglected for the purpose of the 
RG flow, and they do not change the critical properties 
of the model [30]. 

A qualitatively new decimation step, which goes be- 
yond Ref. [30] occurs when the strong bond Ji connects 
two doublet sites. In this case the two sites form a unique 
non-degenerate ground state: 

= 0, rii+i = 1) + \ni ^ 1, Tii+i = 0) 

m,r+i) = , (5) 

which has energy — Jj/2. The second order perturbation 
theory leads to an effective hopping between sites i — 1 
and i + 2: 

Ji-l,i+l = Ji-lJi+l/ Ji = Ji-lJi+l/^- (6) 
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Since each doublet site can be thought of as a spin- 1/2 de- 
gree of freedom, the elimination of consists of the for- 
mation of a singlet. Hence we recover the Ma-Dasgupta 
RG transformation [32, 33]. Note that formally Eqs. (6) 
and (2) are identical. 

C. Flow equations 

Next, we describe the flow equations implied by 
the above decimation steps. As in Ref. [33, 34], we 
parametrize the cutoff energy scale with the variable 
r = log(rio/i^), where fig is the initial cutoff. Also, we 
define the dimensionless couplings Q ^ ^/Ui — 1 which 
are characterized by probability distributions fq{C, T) for 
integer-sites, and /p(C, T) for half-sites. In principle these 
distributions can be different, but one can show that their 
difference is irrelevant in the RG sense, and therefore does 
not affect any of our conclusions. For simplicity, we as- 
sume from the beginning that fq = fp = f. We also 
define Pi = \og{^l/Ji) as the logarithmic bond variable, 
with distribution g(/3,r). Note that by construction f3i 
and Ci have nonzero probability distribution in the inter- 
val [0, oo). 

Renormalization group steps gradually decrease the 
number of remaining sites in the chain (A^(r)). Thus 



decimation of the integer site with large charging gap U 
reduces N by one while a similar decimation of a half- 
integer site simply converts it to the doublet site. Also, 
decimation of a strong link and joining two sites into a 
cluster reduces the number of active sites by one unless 
the link connects two doublets. In the latter case the 
number of remaining sites is reduced by two. Thus the 
flow of N is given by 

dN{r) = - [.go(r)(i + ,s2) + qfo{T)] N{r)dr, (7) 

where /o(r) ee /(0,r) and go(r) = ff(0,r). From Eq. (7) 
and the above RG conditions, we obtain the flows of the 
fractions p, q, and s: 

ds 

^ =-ffos(l-s2) +/o(p + gs) 

^ = -50 [P (1 - s') - 2g (1 - q)] foip -pq) (8) 
^ = -go - 1 + 3g - 2g2 _ qs^] - f„(q - q^). 

It is easy to check that ds/dT + dp/dT + dq/dT = 0, 
provided that p + q -\- s = 1 . 

The RG conditions also lead to master equations for 
the distributions: 



dV 

dgjp) ^ dg{P) 
dV dl3 



dCidC2/(Ci)/(C2)'^(Ci + C2 + 1 - C) - /(C)5o(l + /(C)(/o + 1) 



+ is^9o + qfo) / / dMP2g{f3i)9{l32)S{f3i+f32-f3)+g{f3)go{l-s^)^qg{(3)fo, 



(9) 



Even though these equations look quite complicated, 
the meaning of the each term is straightforward. For ex- 
ample, the second term in the first of these equations 
corresponds to renormalization of the capacitance of the 
cluster following the decimation of the link. The multi- 
plier 1 — s reflects the fact that the renormalization takes 
place only if the link does not connect two singlets. 

The equations (8) and (9) can be significantly simpli- 
fied noting that p = q = (1 — s)/2 is their solution for ar- 
bitrary functions /o(r) and go{T). It is easy to check that 
p ~ q is in fact an attractive solution. Indeed Eqs. (8) 
give: 

^ = 'foil -q)-9o (2(1 -q) + s{l s)) . 

q — p di 

(10) 

Unless q ^ 1, the RHS of Eq. (10) is always negative, 
which means that the line p = g is an attractor. Physi- 
cally one can understand this result as follows: The inte- 
ger versus half integer filling of a cluster is determined by 
the parity of the total number of half-integer decimated 



sites. As clusters grow in size under RG due to coarse- 
graining, the number of such sites becomes large, and 
thus even and odd parities occur with the same proba- 
bility. The case of commensurate disorder, q = 1, which 
was analyzed in Ref. [30] is an exception, since it corre- 
sponds to the strictly zero fraction of half-sites where the 
clusters always remain even. 

The other important observation is that in the weakly 
interacting regime fo^l one can use a simple exponen- 
tial ansatz to solve Eqs. (9): 

As we will see below (sec also Ref. [30]), the univer- 
sal properties of the superfluid-insulator transition are 
determined by the noninteracting fixed point with van- 
ishing /o, where the ansatz is well justified. According 
to our numerical simulations, these exponential scaling 
functions are attractors of the flow equations, and they 
describe very well the distribution of C and /3 even when 
/o is not very small (see discussion below). Substituting 
the ansatz (11) and p = g = (1 — s)/2 into Eqs. (8) and 
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(9) wc find: 



dfo 

dr 

dgo 
dT 
ds 

dr 



^ /o[l-5o(l-s)(l + /o)], 



-f[(l-^)/o + 2.^50] 
^(1-. 2) -50^(1-.). 



(12) 



This system has two fixed points for s: s ~ and s = 1. 
Tlic first fixed point s = 0, p = (7=l/2aswe will see 
below describes the superfiuid phase, while the second 
one: s = l, p = q = corresponds to the random-singlet 
glass insulator. 



D. SF fixed point 

Let us first address the superfiuid fixed point - s = 
(no doublet sites). Note that from the last equation of 
the system (12), this fixed point is stable only when /o 
is small. Since /q fiows either to zero or to infinity, the 
fixed point s = can is stable when /o 0. Then the 
linearized fiow equations to first order in s and /q reduce 
to: 



dT 
dgo 
dT 



-1^ = /o(l -50) 



-2./0. 



(13) 
(14) 



Remarkably, apart from a factor of 1/2 in the second 
equation, these are the same flow equations as obtained 
near the SF fixed point of the random Bose-Hubbard 
model with no half-sites (i.e. (7 = 1). as in Eq. (6) of 
Rcf. [30]. The extra factor of 1/2 can be absorbed in a 
redefinition of /q. This factor appears because half of the 
sites in the lattice have half-integer filling and hence their 
interactions are ineffective in suppressing the superfluid- 
ity and thus do not renormalize go which is related to SF 
stiffness (see Sec. V). 

Equations (13) and (14) are easily solved to give 

/o(r) = A + 2,9o(r) - 21n,9o(r) - 2 = ^ + [1 - .9o(r)]' . 

(15) 

where A is a tuning parameter that controls the flow. 
When A < 0, the flows terminate on the non- interacting 
fixed line fo ~ 0, go > ^■ 

The point A ^ lies on the critical manifold, which 
terminates in a K-T fixed point (note that if we use 
V = y/Jo we obtain the standard Kosterlitz-Thouless fiow 
equations) at 50 = 1, /o = 0. 

At the critical manifold we can solve exatly for the fiow 
of go and fo- Using the parametrization go = I + Sg. fo 
and Sg are determined by the differential equations: 



dSg _ 

or — 

dfo _ 

dr 



4/o 



-foi 



(16) 



By dividing the two equations we have: 
d{Sgf = dfo 

and: 

dSg 
dT 

This gives: 



2 ^ 



5g = T + w 



fo ~ w- 



(17) 
(18) 
(19) 



Note that the solution of the flow equations in the case 
of only integer sites (q=l) yields the same results, only 
with fo ^ 2/r2. 

When yl > 0, the parameters fo and 50 flow past the 
fixed point, where fo begins to increase. The increase of 
fo entails a fiow away from the s = fixed point, towards 
the s = 1 fixed point. In Fig. 1 we show the examples of 
fiows of s near the critical point both in the superfiuid 
and in the insulating regimes. 




50 100 150 200 



FIG. 1: RG flow of the fraction of the singlet sites s according 
to Eqs. (12) near tlie critical point. The solid line corresponds 
to the insulating regime, while the dashed line does to the su- 
perfluid phase. In both cases we assumed there are no singlet 
sites at the onset, i.e. s{T = 0) = 0. 

We now want to stress a very important point. It ap- 
pears that in the superfiuid regime the system fiows to 
the classical fixed line, where there is no charging term. 
However, as we hinted above (see Sec. IIIB), this state- 
ment should be understood with special care. Each time 
the RG scheme merges two sites into a single superfiuid 
cluster it neglects the harmonic Josephson plasmon be- 
tween the two sites, although its energy is well below the 
RG cutoff at this step. These internal excitations do not 
influence the progression of the RG fiow. However, they 
become the elementary phonon excitations of the single 
superfiuid cluster that evolves to be the RG fixed point. 
Therefore the fact that the charging term becomes irrel- 
evant simply means that one should keep only these har- 
monic phonons. In other words, one can ignore vortices 
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or phase-slips, which destroy the superfluid phase if they 
prohferate. Indeed there is a direct connection between 
onsite interactions and phase sHps. As we show next, at 
strong disorder, such sites are responsible for renormal- 
ization of the superfluid stiffness ps playing the role of 
phase slips. The fact that /o flows to zero implies that 
such events renornializing pa become unimportant and 
one can use a nonintcracting quadratic description. This 
issue will be discussed thoroughly in Sec. V. Technically 
the fact that interactions are irrelevant in our description 
comes from the fact that we are working in the grand- 
canonical ensemble. While in the insulating regime there 
is no difference between excitation energies in canonical 
and grand-canonical ensembles, in the superfluid regime 
there is a significant difference. Thus in canonical ensem- 
ble the lowest energy excitation corresponds to a phase 
twist or phonon while in the grand-canonical ensemble 
the lowest energy corresponds to the addition of an ex- 
tra particle, which costs much less energy than the phase 
twist. So the fact that in our scheme the interactions are 
irrelevant in the SF phase should be understood only in 
this grand-canonical sense. 



E. Insulating fixed point - random-singlet glass 

The s = 1 fixed point corresponds to the insulating 
phase. Indeed one can check that in this limit 50 (r) 
and ./o(r) — ^ 00. However, this insulator is not the 
Mott-glass that describes the case of integer only sites 
considered in Ref. [30]. At s = 1, all the sites remaining 
in the system are doubly degenerate. These sites can be 
thought of as spin-1/2 degrees of freedom, with | f) = 
|n+ 1/2) and | [) = |n- 1/2). Without hopping, the 
ground state obviously has a huge degeneracy. However, 
the residual hopping lifts this degeneracy. In the spin 
language, the hopping terms correspond to usual xy spin- 
spin interaction: 



This system implies that /o diverges as /o 



This 



Ji cos ( 



(20) 



Thus we arrive at a spin-1/2 system with random xy 
couplings. The ground state of this system is known to 
be the random-singlet phase [33] . A strong bond between 
sites i and i + \ delocalizes a boson between two sites, 
and creates a cluster that has a charge gap J^: [i/'i, i+i) ~ 
|li)|Oi-)-i) -|- |Oi)|li+i). Quantum fluctuations produce an 
effective coupling between sites i~l and i-\-2 as in Eq. (2). 
The typical length over which the singlets form, i, scales 
as r^. Alternatively, one can say that the gaps of each 
singlet-cluster is A = r2oexp(— V^). 

The flow equations can be linearized near s = 1, and 
we obtain: 



dT 
dga 
dT 
ds 
dT ' 



fo 



-9l 



(21) 



divergence implies that interaction in the remaining non- 
singlet sites is narrowly distributed near the maximum 
energy scale Q,. Following from that scaling of /o, s con- 
verges to 1 extremely fast: s = 1 — (5exp(— e'"). And 
finally go, which corresponds to the average of In(Jj), 
follows the random singlet scaling, and flows slowly to 
zero as: 50 = 

The random-singlet glass is an insulator with the su- 
perfluid stiffness of a chain of length L scaling as [44] 



Ps 



(22) 



with C being a nonuniversal constant. This behavior of 
Ps immediately follows from Eqs. (7) and (21), see also 
Ref. [33]. At the same time this insulator is gapless, with 
the gap also decaying exponentially with y/h but with the 
coefficient C/2. Unlike the Mott glass phase or a Bose 
glass phase, which we will discuss below, the random sin- 
glet insulator is characterized by a diverging density of 
states at zero energy and hence by a divergent compress- 
ibility (k) and superfluid susceptibility Xs- The former 
K = dn/dp, in the spin language is the response to a field 
a^5p. Similarly. Xs defined as the response to the per- 
turbation (5(a + a^), in the spin language, corresponds to 
the perturbation a^6ip. Since the random singlet ground 
state has SU(2) symmetry, the two responses have the 
same form, and diverge as: 



1 



Xs 



Sp\og%no/6p) 



A 



1 



5Alog''(17o/(5A) 



(23) 



(1 - s)fo 



As the slow decay of go with F suggests, indeed the 
random-singlet glass has more superfluid features than 
the Mott-glass. Both have a vanishing gap, but the Mott- 
Glass has vanishing compressibility, and its supcrfluid- 
susceptibility is only finite. 



F. Numerical RSRG for the p-h symmetric state 

In order to corroborate the analytical results for the 
RSRG, we carried out the RG flow numerically with- 
out any simplifying assumptions. The numerics, by and 
large, backs the analytical results. In Fig. 2 several flow 
traces are given in the /o vs. go parameter space, and in 
the p q plane. The initial distributions used consist of 
box distributions in the range 0.2 + 6 < U < 1.2 + S for 
the charging energy, and 0.2 — 5 < J < 1.2 — 6 for the 
nearest-neighbor tunneling. 



IV. PHASE DIAGRAM OF THE B-H MODEL 
WITH GENERIC DISORDER 

When considering experimentally realizable models, 
we must also consider randomness in the chemical po- 
tential, or random offset charges. In particular n could 
have values anywhere between —1/2 and 1/2. Typically, 
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a. 



2.0 



b. 

P 

(half- 
sites) 



0.1 o.s 0.8 1.0 (integer- 
sites) 



FIG. 2: RG flows of various realizations of disorder in the 
(a) /o - go plane, and (b) p-q plane, for chains with ini- 
tial distributions characterized by: {S,po) — (—0.1,0.05) 
(circles), (-0.05,0.04) (squares), (-0.02,0.04) (triangles), 
(0,0.04) (stars). Here po is the initial fraction of the half 
integer p-sites. 



this type of disorder is very relevant. Indeed if we join 
two sites 1 and 2 together into a cluster then the new 
value of ni2 becomes a sum of ni and n2 modulo one (so 
that the result also belongs to the (-1/2, 1/2] interval). 

To address this problem it is worthwhile review a few 
important insights gained from our analysis of the par- 
ticle hole symmetric model. In that case the low energy 
behavior was dominated by the line p = q where the 
number of half-sites with n = 1/2 was the same as the 
number of integer-sites, with n = 0. At the same time 
we saw that the universal properties of the SF-INS tran- 
sition with p = q and p = 0, q = 1 were identical up to 
a factor of one half in Eq. (14), which is absent in the 
integer filling model. 

We thus can anticipate that the fixed point governing 
the SF-insulator transition has a uniform distribution of 
n when we remove the particlc-holc symmetry restric- 
tions. But also, in analogy with the half integer case, 
we can expect that the fixed point describing the SF-INS 
transition remains intact. On the other hand, again hav- 
ing the p-h symmetric case in mind, we expect that the 
distribution of n at the critical point strongly affects the 



properties of the insulating phase. As it turns out, the di- 
agonal disorder plays the role of a 'dangerously irrelevant 
variable' (as s is in the analysis above - irrelevant in the 
SF side of the transition but strongly relevant in the in- 
sulator side) . The diagonal disorder does not change the 
nature of the critical and crossover behavior, but it de- 
termines to which insulating phase the system will flow. 
In Ref. [30] the p-h symmetric model with only integer 
fillings had a Mott-glass insulating phase. By allowing 
also sites with charging degeneracy {n = 1/2) but pre- 
serving the p-h symmetry, the system in its insulating 
phase is a random-singlet glass. When we remove the 
p-h symmetry, we expect that the insulator becomes a 
Bose-glass: gapless, compressible state with a diverging 
susceptibility to SF fiuctuations. 

Before going into more detailed analysis, which con- 
firms the above assertions, we would like to comment 
on the similarities with the perturbative RG approach 
of Giamarchi and Schulz [25]. In particular they derived 
the following fiow equations near the transition between 
superfiuid and localized phases: 



dr 2 



dr 



1 

2^ 



(24) 
(25) 



where K is the Luttingcr parameter, ^/pJk and V is pro- 
portional to the variance of the disorder in the chemical 
potential: ii{x)fi{x') cx 'D5{x — x'). We point out that 
in Ref. [25] the Eqs. (24-25) are written in terms of the 
inverse Luttinger parameter being the inverse of the 
common convention which we use (see e. g. Ref. [37]). 
Note that there is a direct analogy between Eqs. (13), 
(14) and Eqs. (24) and (25) if one identifies K with qq 
and T) with /q. 

Interestingly in Ref. [27] it was argued that the SF- 
INS transition described by Eqs. (24) and (25) does not 
belong to the KT universality class because of the first 
power of T) appearing in the second of these equations, 
as opposed to the second power in the conventional case. 
This difference according to the authors lead in particu- 
lar to the unconventional scaling of the correlation length 
with K — Kc near the critical point. However, this must 
be a misstatement since, as we argued earlier, the sub- 
stitution V = brings the flow equations to the con- 
ventional KT form. This change of variables should not 
affect the scaling. Moreover the flow equations in terms 
of K and a are more natural because a has dimensions of 
the external potential and thus it (not V) is analogous to 
the strength of the periodic potential, which drives the 
transition in a nondisordercd case. 

The similarity between the perturbative analysis of 
Ref. [25] and the one presented here goes even further. A 
simple scaling argument shows that disorder in the chem- 
ical potential is strongly relevant in both approaches. 
However in the language of Ref. [25] the strongly relevant 
part of the disorder corresponds to the forward scatter- 
ing, which can be reabsorbed into the canonical smooth 
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fluctuations of the density. It is the backward scattering 
or phase-slips, which determine the fate of the superfluid 
phase. By analogy with our approach we can argue that 
even the smooth part of the disorder potential should 
become strongly relevant in the insulating regime. Thus 
it should play the role of a dangerously irrelevant term 
just as the distribution of n docs in our approach. Unfor- 
tunately the pcrtubativc RG approach becomes uncon- 
trolled in the insulating regime and this postulate cannot 
be reliably verified. 



A. RG scheme for the generic disorder B-H model 

We probe the observsations above by extending our 
RSRG analysis to treat arbitrary disorder: C/^, Ji, and rii 
all random in the model in Eq. (1). 

First let us analyze the charging term while ignoring 
the hopping. Each site has a charge gap given by: 



1. 



:Udl~2\n, 



(26) 



where -1/2 < n< 1/2. 

As before, we treat the model iteratively, but this time, 
in each step of the RG we find the largest energy scale: 



= max{ Ji, Ai} 



(27) 



and eliminate it. If it is a gap, A^, then the site i freezes 
into its lowest energy charging state. Quantum fluctua- 



tions induce an effective hopping between sites i - 
i - 1: 



J, 



Ji-lJi 



1 



n (l + 2|n,|)- 



1 and 



(28) 



The last multiplier in this expression is a non-universal 
prefactor, which varies between 1 and 1/2. This prefactor 
does not affect any universal features of the transition 
and we can safely set it to unity. On the other hand, if 
Ji is the largest energy scale, then the sites i and i + 1 
form a SF cluster with effective charging energy given by 
Eq. (3), and with a filling offset: 



H, i+1 



"2 



(29) 



where the last equality is defined modulo adding or sub- 
tracting one, so the the result always belong to the inter- 
val: ni,,+i e (-1/2,1/2]. 



B. Generic case flow equations 

The ensuing flow can be quantified using flow equa- 
tions for the distribution of logarithmic couplings, 
(3i = logfl/Ji, and the joint distribution F{(,n) = 
f{C,n)x{C - 1 + 2|n|) where Q = Q/U and the Hcavi- 
side step function x enforces the constraint fi/A > 1 or 
equivalently C > 1 — 2|rT|. 

The flow equations are given by: 



- ll + [A + .9o(i - Is)] 9 X g^^^^+p^ + gigofs - A) 



dT dp 
dl 
dT 



df df 

+ .90 / X /|^^^^^^^ + /(I + /i - gofs), 



(30) 
(31) 



r 



0.5 



where 



fi = J dn(l-2|n|)/(l-2|n|,n) (32) 

-0.5 

is the density of sites with a large charging energy, 

fs-jdcj dnxiC - 1 + 2\n\) f X /|^^^^^^^ ^_ . 

(33) 

Physically 1 — /s is the density of strong bonds connecting 
the sites with large onsite interaction, which are close 
to half filling. These sites form a cluster with A > 
and thus have to be eliminated as a spin singlet. In the 
equations above x implies the convolution over /? in g x g 
and over both and n in / x /. 



Although the equations look somewhat obstruse, near 
the critical point they can be solved with the same scal- 
ing ansatz as before. Indeed, since near the transition the 
interactions are negligible, we can safely ignore 6C, which 
is of the order of one, in the convolution (33). Also sim- 
ilarly to the particle-hole symmetric case we can expect 
that near the critical point the distribution of n is uni- 
form and thus f{C,ri) does not depend on n. We then 
use our standard scaling ansazt: 



g(/?) = 5oe 
/(CrT) 



-ffo/3 



1 



-/oC 



/oe" 



'/oC 



(34) 
(35) 



where in the last equality we used /o <C 1. In the same 
approximation of small /q wc find that /i w /o/2. Substi- 
tuting the scaling ansatz into the flow equations (30) and 
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FIG. 3: Flows in the /o vs. go plane. Initially, the interaction 
energy U and the hopping J are uniformally distributed in 
the range 0.2-5 < J < 1.2-5, and 0.2 + S <U < 1.2 + 5; the 
offset charge is also uniformally distributed between — An < 
n < An. The values of 5 and An for the plots shown are: 
S = —0.05, An — 0.12 which is in the SF phase (circles), 
5 = -0.04, An = 0.08 (squares), 5 = -0.03, An = 0.08 (tri- 
angles), S = -0.02, An = 0.08 (stars), and 5 = 0, An = 0.12 
(diamonds) . 



(31) and using /o ^ 1 wc immediately recover that /o 
and go obey Eqs. (13) and (14). The latter automatically 
implies that the SF-IN transition in the case of generic 
disorder belongs to the same universality class as in the 
particle-hole symmetric case. 

We confirm these findings performing numerical anal- 
ysis of the full RG equations (30) and (31). We find a 
clear signature of the K-T transition that is even more 
pronounced than before. In Fig. 3 the flows in the f-g 
parameter space are shown for three different initial con- 
ditions. 



C. Nature of the Bose-glass 

We now turn to the analysis of the insulating phase at 
generic disorder, i.e. of the Bose glass. Even though the 
simple exponential form of fiCn) and g{(3) does not give 
the exact solution to the flow equations, as we deduce 
from numerical analysis, it gives a very good approxima- 
tion to the true distributions. For large values of /o the 
equation (35) simplifies to f{C,,ri) « /o exp(— /oC). It is 
straightforward to check that under these conditions we 
have /i ~ 1 and fs ~ 4/o exp(— /o/2). Then the flow 
elquation for go becomes very simple: 



dm 

dV 



-go- 



(36) 




FIG. 4: A semi-log plot of go vs. F in the Bose glass phase. 
The two samples are the same as the star and diamond curves 
in Fig. 3 with: S = -0.02, An = 0.08 (stars), and S = 
0, An — 0.12 (diamonds). The dark lines are guide to the 
eye, and have a slope of —1 in the plot. As can be seen, the 
late stage of the flow of go fits go ~ e""" very well. 



and the random singlet insulator with g'^ sa —go- Phys- 
ically the parameter go characterizes the strength of the 
hoppings remaining in the system. As we argued before 
slow 1 /F decay of go in the random singlet phase resulted 
in the divergent density of states at zero energy and as 
a result in a divergent compressibility. On the other 
hand in the Mott glass go was vanishing very rapidly 
go ~ exp— (exp(F)) and thus the corresponding Mott 
glass had a vanishing density of states at zero energy 
and vanishing compressibility. In the Bose glass phase 
we have exactly intermediate behavior: go ^ exp(— F). 
As we will see shortly this scaling implies finite density 
of states at low energies and thus a finite compressibil- 
ity. In Fig. 4 wc plot the flows of the parameter go as 
a function of F obtained from numerical solution of the 
RG equations for two different insulating samples with 
generic disorder. As can be seen from the latter figure, 
the flow of go vs. F is consistent with go ~ exp(— F). 

Another interesting conclusion is coming from the fact 
that /i « 1 is independent on fo (and hence dgo/dV « 
—go)- We remind that /idF gives the probability of dec- 
imating a site when we change the cutoff scale from to 
il(l — dr). It is remarkable that in the Bose glass phase 
this probability is independent of fo- On the contrary, 
the probability of eliminating the link is proportional to 
go and thus vanishes at long F. Thus wc come to the 
conclusion that the number of the sites remaining in the 
system TV scales exactly as the cutoff energy scale: 



N 



(37) 



Such flow indicates that the Bose glass phase is indeed 
intermediate between the Mott Glass where w —fo9o 



which indicates uniform density of localized states in the 
insulating state. As we will see shortly the parameter 
K plays the role of the compressibility. It is interesting 
to note that k discontinuously changes across the phase 
transition. Indeed very close to the transition the number 
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of sites remaining in the system is given by 
N{T) ^ iVoexp(-.gor) ^ exp(-r). 



(38) 



This behavior is correct for length scales shorter than 
the correlation length ^ ~ exp(l/-\/^), where A is the 
tuning parameter appearing in Eq. (15), A corre- 
sponding to the critical point. After that we should use 
the flow equations valid in the insulating regime where 
N ~ exp(— /iF) exp(— r). So we see that the scal- 
ing N ^ Nq exp(— r) works very well. We thus conclude 
that the ratio iV(r)/il(r) goes to a constant, which is 
independent of A. 



1. Compressibility and SF susceptibility 

The easiest way to see that k is indeed the compress- 
ibility is to map the renormalized array of clusters into a 
spin-1/2 chain. Since deep in the insulating phase the 
displacement \ri\ is close to 1/2, the local interaction 
strengths, t/^, is mostly quite large, and obeys: Ui > ft. 
This implies that we could retain only the two lowest 
charing states, which we the spin up and spin down states 
of an effective spin 1/2 degree of freedom: 



1 



(39) 



In this picture the gap = U{1 — 2\n\) plays the role 
of the external magnetic field along the z axis /if. The 
hopping Ji is in turn maps to the xy coupling between 
the neighboring spins. 

Let us first determine the distribution function of h^: 
H{h^) assuming that f{C,n) is given by Eq. (35) with 
/o » 1. 

1/2 oo 

Hm^fS jdn j dCe~f^'^5[\h^-\-^{l-2\n\) 



2Vt 



'1/2 l-2|rr| 



(40) 



which is just a uniform distribution. 

The compressibility of the insulating phase is given 
by the z-field susceptibility of the spin chain. The lat- 
ter is easily shown to be twice the probability density of 
H{h'^ = 0). Thus the compressibility is: 



dn 



N 



2NH{h^ = 0) 



N 



(41) 



where h^^^ is the infinitesimal external magnetic field 
along the z axis. Indeed, this is the result advertised 
inEq. (37). 

The superfluid susceptibilty Xs is obtained as the re- 
sponse of the spin chain on a small magnetic field in x 
direction. Note that in the Bose glass phase the coupling 
between different sites is vanishingly small, and thus Xs 



can be derived by considering an isolated site, which is 
described by a spin 1/2 Hamiltonian: 

= h's~ + h-Js^. (42) 

A straightforward calculation yields that the average 
magnetization along the x axis is: 



1 h'^ 

_ '^cxt 

Thus the susceptibility is: 

f d/i^ 1 N n 

Xs = N ^ log 



(43) 



2n \h'' 



2^ ^\hU 



K 17 



(44) 

Obviously Xs diverges as h^^^ 0. Thus we find that 
the Bose glass phase is characterized by divergent ps and 
finite k in agreement with Ref. [27]. 



V. COMPRESSIBILITY, STIFFNESS, AND 
LUTTINGER PARAMETER IN THE 
SUPERFLUID PHASE AND AT CRITICALITY 

Let us now focus on the properties of the superfluid 
phase we find. The superfluid phase is associated with 
the formation of a superfluid cluster that spans the chain. 
The cluster consists of all the original bare sites that 
were not decimated due to their charging energies. These 
surviving sites obey the Hamiltonian: 



njSF 

^eff 



E 



\2 



(45) 



where ipi and n,; are the phase and number operators of 
the surviving cites, and Ui are the charging energies of 
each bare site. The Ji are harmonic couplings between 
the surviving sites, which are the result of the decimation 
of a strong bond (marked with a tilde since they can get 
renormalized by intervening charge-blockaded sites) as 
we now explain. When deriving the RG equations, we 
eliminated the strongest bonds iteratively, by setting the 
sites they connected into phase-coherent clusters, which 
implies replacing the strongest Josephson couplings with 
a harmonic coupling: 

1 2 

-Ji cos {ipi+i - Lpi) -Ji {ipi+i - ipi) . (46) 

We then approximated the cluster to be phase-coherent: 



f~ (ft 



V2- 



(47) 



This strong approximation is sufRcient for obtaining the 
flow equations, but we need to allow intra-cluster fluctua- 
tions it in order to discuss the properties of the superfluid 
phase. 
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The stiffness and the compressibihty of the supcrfluid 
phase are given in terms of the parameters Ji and Ui 
in the effective Hamiltonian (45), which describes the 
prohferating superfluid chister. The compressibihty is 
given by 



We now use that: 



1 



LUsF- 



1 



cluster 



-E- 



(48) 



where L is the total length of the chain, 
superfluid stiffness is similarly obtained as 



1 

Ps 



The inverse 



(49) 



Note that the simple expression for the stiffness owes to 
the fact that the fixed point Hamiltonian (45) is har- 
monic. Therefore the stiffness suffers no further renor- 
malization by quantum fluctuations and it is the same as 
in the classical model (see [30] ) . 

We will now proceed to calculate the average compress- 
ibility , stiffness, and Luttinger parameter of the super- 
fluid K = IXjKps . 



A. Differential equation for the inverse charging 
energy 

The compressibility given by Eq. (48) can be calcu- 
lated in a rather straight forward way within the RG 
scheme outlaid in the previous sections. The variable C 
in the RG scheme is specifically designed to keep track 
of the cluster compressibilities. We recall the RG flow 
equation (9) for the distribution function /(C) 



d/(C) 



(1 ■ 

+90 J dCidC2S{C - Ci - (-2 ' I) 
/(Ci)/(C2) + ./(C)(./o + l-5o). 



(50) 



The solution to this equation will allow us to compute the 
compressibility from the average value of C as k = C/^- 
To obtain a differential equation directly for the av- 
erage compressibility (inverse-charging energy) we move 
to the Laplace-transformed representation: F{ri) = 



e ''''/(C), which obeys: 



^ = -fo + F{^){rj-l)-v^ 
+.goF(r;)2e-" + F(77)(/o + l-go). 



(51) 



dF{r]) 



(52) 



to obtain: 



dr 



C(/o + .go - 1) - 1 + 50 



(53) 



The inverse charging energy is given by Eq. (48), in 
which an extra factor of appears. Adding it on we 
obtain: 



dr 

where 7 = C/^- 



7(/o + 5o)-(l-5o)/f^ 



(54) 



B. Flow equation for the stiffness 

Calculation of the stiffness requires a slight extension 
of the RG scheme. The method described thus far did 
not include a cluster variable which stores the internal 
stiffness. In other words the RG scheme does not keep 
track of the internal sum over 1/Ji (49) within the pro- 
liferating clusters. 

Fortunately, such a variable can easily be easily in- 
cluded by extending the cluster distribution function /(C) 
to a joint distribution /(C,x)) where 



^ 7 

i^cluster * 



(55) 



is a variable designed to keep track of the superfluid stiff- 
ness of the clusters. Each time two clusters are joined in 
the RG flow by a large bond J — fl, the variable x of the 
joined cluster is given by: 



X12 = Xi + 



Xi+X2 + -^e^- (56) 

"0 



The flow equation for /(C, x) is a- straight forward ex- 
tension of Eq. (50): 



dfiC-X) _ 



9C 



= (1 + 0- 

+50 / dCidC2SiC - Ci - C2 - 1) / dxidx25{x - Xi - X2 - ^o'e^)fiCi,Xi)f{C2, X2) 

+/(C,x)(/o + l-5o). 



(57) 
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This is a rather compHcatcd equation for the joint dis- 
tribution of cluster stiffness and charging energy. How- 
ever it can be greatly simplified if we are interested 
only in the average of the stiffness. The latter can 
be calculated by integrating Eq. (57) with respect to 
C, and taking its Laplace transform with respect to x- 

oo oo 

^(A) ^JdCj dxe-^^fiC, X)- This yields: 





dS{\) 

dr 



= -/(0,A)-5(A) 
W-pr + 5(A)(l + /o 



Where: 

/(0,A) =/dxe-^^/(C = 0,x) 
Again using the fact that: 



(58) 



.90 j 



we obtain: 



dr = ^^-^^ 



X 



go) + ga/Q 



dS{X) 




dX 





9/(0, A) 



dX 



(59) 



(60) 



where we used = fig exp(— F) We note that the only 
difference between Eq. (54) and Eq. (60) is in the sub- 
leading term, 50/^ above, and {go — l)/ri in Eq. (54). 
There is also the last term in Eq. (60), which should be 
negligible and negative. 



C. Differential equation for the length of a 
superfiuid cluster 

The differential equation for the typical length of the 
clusters [30] is given by 



.90) 



(61) 



It is interesting to note that this equation is the same, 
at the leading order, as the equations, derived above, for 
the sums of inverse charging energies ^ C s-^d the sum of 
inverse Josephson couplings x within a cluster. 

When calculating the stiffness and compressibility us- 
ing the flow equations for a particular cluster, as illus- 
trated above, we need to renormalize until the size of a 
SF cluster is that of the entire chain: 



£r = L. 

Therefore the compressibility: 

1 



1 

K = — - 



L UsF~cluster 

and the inverse stiffness: 
1 

Ps 



-E- 



(62) 
(63) 
(64) 



D. Compressibility and stiffness at the critical 
point 

Let us assume that we start sufficiently close to the 
critical point, so that the distributions /(C, F) and g{(3, F) 
already converged to the universal forms characterized by 
/o(r) and go(r) (see Eqs. (34) and (35)). In Sec. (HID) 
we found the explicit flow of these functions at criticality 
.90 - 1 « 2/F + 2/F2 and fo « 4/F2 (that is, the flow on 
the separatrix). These flows start at some initial value 
Fq which characterizes the bare disorder distributions of 
the microscopic system. The larger is Fq the wider is the 
disorder distribution in J^. 

Combining this with flow equation for C (53) we find 



dT M F r2 



2 

F' 



This has the solution for large F: 

C«Cc-F2exp[-6/F] 



(65) 



with being a constant, which we obtain from initial 
conditions. We know that for sufficiently large Fq, Co ^ 
= iF^. This implies: 

Jo 



1 



C=-F^exp[-6/F] 

and = 1/4. 

Similarly, for x we get: 

X « ^^r2 exp[-6/F] {C^ - exp[6/F]) , 



(66) 



(67) 



where the second term in the brackets comes from the go 
term in the equation for x- The constant can also be 
obtained from boundary conditions. At the onset Xo = 
since we start with bare sites, and only after some RG 
we get the sum of 1/J to grow. This implies: 

(68) 

By the same token, the solution of Eq. (61) at criti- 
cality is: 



Cx = exp[6/Fo]. 



^ 



(69) 



where £0 is of order 1. 

Putting our results in the definitions of ps and k, we 
obtain the compressibility of the critical system sending 
F 00: 



and the inverse stiffness: 

1 = 1 = - 

Ps £r 



F2 
^ 



^ 

Afloio 



■(exp[6/Fo] - 1) 



(70) 



(71) 



always tend to a number as F — !■ 00. 



The energy scale for both the stiffness and the compress- 
ibility is given by Qq, the initial energy scale of the prob- 
lem. Both also tend to constants along the critical flow 
line. 
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E. Luttinger parameter at criticality 

By multiplying Eqs. (70) and (71) we obtain the lut- 
tinger parameter of the SF cluster: 

K' = ^np,= ^-^ I (72) 
2 exp[6/roJ - 1 

Indeed we find that it is a constant along flows on the 
critical manifold, which is independent of the initial en- 
ergy scale fio- On the other hand this result is clearly 
not universal, since it depends on Fg. 

As mentioned above, Fq parameterizes the strength of 
the bare bond disorder distribution. For a given system 
on the critical manifold, the larger is Fq the broader is 
the system's initial distribution of both J and 1/U . We 
can therefore interpret Eq. (72) as stating that at strong 
disorder, the Luttinger parameter required to stabilize 
a superfluid phase depends on the disorder strength. A 
larger Luttinger parameter is needed the more disordered 
is the system. This statement is clearly different from 
the situation at weak disorder, for which Giamarchi and 
schulz had predicted a transition at a universal value of 
the Luttinger parameter [24, 25]. We shall comment on 
the relation between these two limits in the discussion 
below. 



VI. DISCUSSION 

Using the real-space RG approach, we obtain a consis- 
tent picture both of the possbile insulating phases of the 
random-Bose-Hubbard model, but also of the transition 
from the superfluid to them. The seminal work of Gia- 
marchi and Schulz (GS)[24, 25] obtained a description of 
what seems to be the same transition in terms of a pcr- 
turbative RG in weak-randomness - the opposite limit to 
our starting point. We now ask: how do these two sce- 
narios, or descriptions, correspond to each other? Now 
that we obtained our result for the Luttinger parameter 
at criticality, Eq. (72), we can address this question. 

One of the central results of Ref. [25] is the universality 
of the Luttinger parameter at the transition: 

Kf'^ = I (73) 

Since GS considered the anomalous dimension of what is 
essentially a phase-slip operator, the universality of K at 
the transition was deduced from the fact that when K > 
3/2, phase slips are irrelevant. Since in weak randomness 
phase slips are clearly the most relevant operators, the 
vanishing of their scaling dimension implies criticality. 
Also, the generality of the GS approach, and the self- 
averaging of the SF phase [45] implies that phase slips 
turn relevant when K = 3/2 even for strong disorder. 

At strong randomness, however, we find that a differ- 
ent type of disturbance of the superfluid phase can dis- 
order it. In the real-space RG analysis grains with large 



charging energies are decimated, implying that a whole 
grain becomes isolated from the rest of the chain. This 
process is equivalent to a phase-slip dipole happening 
around the grain. Phase-slip dipoles consist of a phase 
slip and an anti-phase slip happening simultaneously at 
neighboring positions in the chain. In the week coupling 
limit, these dipoles are not enough to degrade the super- 
fluidity, since they do not produce a voltage drop. But 
when the disorder is strong, the dipoles, or equivalently, 
the blockaded insulating sites, suppress tunneling across 
the lattice, as we find from our analysis. 

For sufficiently strong disorder, the Luttinger param- 
eter at which blockaded sites destroy superfluidity, i.e., 
the critical Luttinger parameter, is given by Eq. (72): 

^ = V^ cxp[6/Fo]-l - ^''^ 

For the p-h symmetric case considered in Sec. Ill, and 
K = 7r-y2/[exp(4/Fo) — 1] for the commensurate case, 
with rij = 0. As explained below Eq. (72), Fq is a 
measure of the initial disorder of the system. Thus, K 
grows monotonically with the disorder, and exceeds the 
universal GS value of K = 3/2 at intermediate values 
of Fq ~ 3. This implies that the transition we find 
takes over the universal GS transition at a finite disorder: 
since we find that the breakdown of superfluidity occurs 
at > 3/2, the transition happens well into the region 
where single phase-slips are irrelevant, and thus they do 
not modify the critical properties of the model, and can 
be safely ignored. This also justifies our procedure of 
SF cluster formation as outlined in Eq. (46). It is in- 
teresting to note that Fq ^ 3 corresponds to a charging 
distribution which is peaked at about 4/Fq51o ~ l/2fioj 
as obtained by plugging Fg into Eqs. (11). 

Our conclusion is that at finite randomness the critical 
fixed point of the RSRG takes over (Fig. 5). When this 
happens, universality of the Luttinger parameter at the 
transition is lost. Since the transition we are describing 
is still a Kosterlitz-Thouless type transition, many prop- 
erties of the weak-randomness transition, and strong ran- 
domness transition are shared. One can argue that the 
Luttinger-parameter universality lost at strong disorder 
morphs into a different universlity — that of the expo- 
nent with which the distribution of J vanishes at small 
energies, which is 50 ^ 1 ^ 0+ at criticality. 

An outstanding question is how the weak-randomness 
phase-slip driven transition changes into the transition 
we find at strong disorder. One possiblity is that the 
two scenario continuously morph into each other. Yet 
another more exciting possiblity is that our analysis is 
equivalent to the calculation of the scaling-dimension of 
an operator different from single phase-slips, and that 
such an operator becomes relevant at sufficiently strong 
disorder at Luttinger parameters K > 3/2. Therefore it 
causes a break down of superfluidity before phase-slips 
become relevant. 

Another important difference between the perturba- 
tive approach of GS and our results is that GS assume 
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weak randomness 
transition 



strong randomness 
transition 



Insulator 



Randomness 

FIG. 5: From all the analyses we carried out it seems that our 
transition does not happen at a universal value of the luttinger 
paerameter, but rather, at a universal value of the power-law 
of the J distribution {go = 1). From the discussion, it seems 
that there are two scenarios for the breakdown of the SF. 
At weak randomness it is the G-S single-vortex proliferation 
mechanism that first destabilizes the SF. In this range, the 
scenario we present would destabilize the SF at lower K than 
single vortex proliferation (dashed grey), and is therefore not 
a true boundary. At larger randomness, our scenario is the 
first to stabilize the SF, as it occurs at larger K's than the 
universal G-S value. A concequence is that the universality 
of the Luttinger parameter at criticality is lost. 



that the diagonal disorder is gaussian and fully charac- 
terized by its variance while the off-diagonal disorder is 
weak and irrelevant. On the contrary, in the strong ran- 
domness approach, we see that transition corresponds to 
a wide power-law distribution of tunneling amplitudes. 
Standard replica methods are not applicable to this type 
of disorder distribution and thus it is not surprising that 
the jump we find in K is different. The appearance of 
broad power law distribution of links in ID is not sur- 
prising. There is always a finite chance of encountering a 
large insulating cluster separating two superfluid regions, 
which effectively blocks the tunneling between superfiu- 
ids. This is a special property of ID systems. In Ref. [30] 
we demonstrated that this is indeed the case for a simple 
toy model. The real space RG just reflects this property 
of ID systems. Thus if disorder is not very small so that 
weak links necessarily occur with finite probability, we 
believe that our scenario of the SF-IN transition to be 
more plausible than GS scenario of weak disorder. How- 
ever, the final resolution of this question is currently be- 
yond reach of both the RSRG and the GS analysis, since 
it is concerned with the intermediate randomness regime. 
Probably this question can be addressed numerically. 



VII. EXPERIMENTAL CONCEQUENCES 

A. Critical current of a finite superfluid chain 

At low energy scales, the real-space renormalization 
group allows detailed knowledge of the superfluid phase. 
Most importantly, the effective low-energy Josephson 
junction coupling distribution is: 



PiJ) = go 



I 



Jl-So ■ 



(75) 



The knowledge of the Josephson distribution function al- 
lows us to make a connection with a rather simply mea- 
surable experimental property: The critical current of a 
chain. 

Unlike the Luttinger parameter, the critical current of 
a bosonic chain in the absence of phase fluctuations is 
controlled by its weakest hopping link. Given a strongly 
disordered bosonic chain in its superfluid phase, we can 
apply the real-space RG until the effective coupling dis- 
tirbutions approach their universal behavior, and in par- 
ituclar the distribution (75) for the Jospehson energies of 
each bond, and with negligible charging effects. 

Let us now calculate the scaling of the critical current 
on the bare length of the system. Given a particular 
disorder distribution, the universal distributions are ob- 
tained once the UV cutoff is ^Iq, and only a fraction 1/i 
of the chain is still active, and the chain is of length L/£. 
The scaling behavior of the weakest Josephson energy ex- 
pectation value, Jmin, is obtained by requiring that the 
probability of having at least one bond with an energy 



L r dj 



90 
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Carrying out the integral we obtain: 

1/90 

L. 



(76) 



(77) 



In the weak disorder regime, where go ^ 1 we see that 
the critical current is almost size independent. While at 
strong disorder near the transition go ^ 1+ the critical 
current scales as the inverse system size. This prediction 
can be directly tested in experiments. Using extreme 
value statistics one can even flnd the whole Gumbel dis- 
tribution of the critical current in the SF regime: 



9o 



90 
min 
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90 



(78) 



B. Resistance at finite temperatures 



By a similar argument, we can guess the finite temper- 
ature behavior of a disordered superfluid chain. First, we 
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make the following simplifying assumptions: if a bond 
strength is J > T, we can neglect its finite temperature 
resistance, but if J < T, a bond will give a finite resis- 
tance r, which is T independent. Furthermore, we ignore, 
for the sake of this discussion, the dependence of r on J. 

Under these simple assumptions, the resistivity p at 
temperature T is given by the density of bonds of 
strength J < T. Therefore: 



dJ 



go 



r{T/nor (79) 



where, as defined above, Qq is the rough energy scale at 
which the chain is exhibiting the universal low energy 
behavior. 

In finite chains, we expect that Eq. (79) would only 
be valid when T S> Jmin- Very crudely, by replacing the 
lower limit of the integral in Eq. (79) by Jmin as given 
by Eq. (77), we obtain for T > Jmin- 



T 



(80) 



VIII. CONCLUSIONS 

In this paper we extend the real-space RG analysis of 
Ref. [30] to the case of non-commensurate chemical po- 
tential. We find that remarkably, the symmetry and de- 
tails of the diagonal disorder are irrelevant for the SF-INS 
transition in a system with only onsite interactions. Nev- 
ertheless, the symmetry of the disorder completely deter- 
mines the type of insulator that the system obtains. The 
superfluid phase will break down at a Kosterlitz-Thouless 
critical point, and will become: (i) a gapless, incompress- 
ible, Mott-glass if the chemical potential is commensurate 
{rij = 0), (ii) a gapless, compressible Bose-glass with di- 
verging superfluid susceptibility if 1/2 < rlj < 1/2 is un- 
restricted, (iii) a gapless random-singlet glass with a di- 
verging compressiblity and superfluid susceptibilty in the 
case of p-h symmetric chemical potential {nj = 0, 1/2). 

An important question about our approach is its 
connection with the seminal work of Giamarchi and 



Schulz [24] , we calculated the properties of the superfluid 
phase using the real-space RG analysis. By considering 
the Luttinger-parameter we showed that at strong 
disordered the SF-INS transition occurs at a finite value 
of K, larger than the universal GS value, and that the 
universality of the Luttinger parameter is replaced with 
a universality of the power-law distribution of effective 
hopping at low energies. The real-space RG approach 
is thus not complementary to the GS approach, but pro- 
vides a description of the SF-INS transition at strong dis- 
order, and allows direct access to the insulating phases, 
where the GS approach fails. 

An interesting direction to pursue in the future is 
the utilization of the RSRG approach for calculation of 
transport-propoerties, and finite-temperature properties 
of the random 1-d Bose-Hubbard chain. Such calcula- 
tions could probably be done by combining our approach 
with that of Motrunich, Damle, and Huse [46]. The pres- 
ence of very large disorder in the insulating phases should 
make such calculations accessible. On the other hand, 
they may prove difficult near the transition due to the 
finite randomness there. 

Another outlying question is that of the correlations 
in the SF phase. Self averaging indicates that the Lut- 
tinger parameter we find in Sec. V also dictates the de- 
cay of correlations in the strongly-disordered superfluid 
phase. This, however, remains to be confirmed in direct 
numerical investigation of a strongly random harmonic 
chain. We should emphasize that due to our method for 
finding the Luttinger parameter, it should be consistent 
with the anomalous dimension of the phase-slip operator 
in the GS theory. 
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